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Abstract 

We consider in this paper a solvable model for the motion of molecular motors. 
Based on the averaging principle, we reduce the problem to a diffusion process on 
a graph. We then calculate the effective speed of transportation of these motors. 
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1 Introduction 

One of the possible ways to model Brownian motors /ratchets is to describe them as 
particles (which model the protein molecules) traveling along a designated track (see [7]). 
At a microscopic scale such a motion is conveniently described as a diffusion process with 
a deterministic drift. On the other hand, the designated track along which the molecule 
is traveling can be viewed as a tubular domain of some random shape. In particular, 
such a domain can have many random "wings" added to it. (See Fig.l. The shaded 
areas represent the "wings".) In this paper we are going to introduce a mathematically 
solvable model of the Brownian motor and discuss some interesting relevant questions 
around this problem. Our model is based on ideas similar to that of [5] and [2j Chapter 
7]- 

The model is as follows. Let h^ix) be a pair of piecewise smooth functions with 
/iq~(x) — Iiq(x) = Iq{x) > 0. Let Do = {( x i z ) '■ x £ M, /iq(x) < z < (x)} be a tubular 
2-d domain of infinite length, i.e. it goes along the whole x-axis. At the discontinuities of 
h^(x), we connect the pieces of the boundary via straight vertical lines. The domain Dq 
models the "main" channel in which the motor is traveling. Let a sequence of "wings" 
Dj (j > 1) be attached to Dq. These wings are attached to Do a t the discontinuities of 
the functions h (x). 
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Consider the union D = Dq {J Dj. An example of such a domain D is shown in 

3=} 

Fig.l, in which one can see four "wings" D\, D2, -D3, -D4. We assume that, after adding 
the "wings", for the domain D, the boundary dD has two smooth pieces: the upper 
boundary and the lower boundary. Let n(x,z) = (ni(x, z), ri2(x, z)) be the inward unit 
normal vector to dD. We make some assumptions on the domain D. 

Assumption 1. The set of points i£l for which there are points (x, z) G dD at 
which the unit normal vector n(x, z) is parallel to the x-axis: n%{x % z) = has no limit 
points in R. Each such point x corresponds to only one point (2, z) G dD for which 
n 2 (x, z) = 0. 

Assumption 2. For every x the cross-section of the region D at level x, i.e., the 
set of all points belonging to D with the first coordinate equal to x, consists of either 
one or two intervals that are its connected components. Also, for some < Iq < Iq < 00 
we have lo < h^{x) — h^{x) = lo(x) < Iq. 

Let us take into account randomness of the domain D. Keeping the above assump- 
tions in mind, we can assume that the functions h (x) and the shape of the wings D]~ 
(k = 1,2,...) are all random. Thus we can view the shape of D as random. We introduce 
a filtration —00 < s < t < 00 as the smallest c-algebra corresponding to the shape 
of D n {(x,z) : x G We introduce stationarity and mixing assumptions. Let us 

consider some 4 £ J 7 ', —00 < s < t < 00. The set A consists of some shapes of the 
domain D C\ {(x,z) : x G [s,i]}. Let 9 r (r £ R) be the operator corresponding to the 
shift along x-direction: 6 r {A) G J^lXr consists of the same shapes as those in A but 
correspond to the domain D n {(x, z) : x G [s + r, t + r]}. 

Assumption 3. (stationarity) We have P(^4) = P(9 r (A)). 

Assumption 4. (mixing) For any A G J- l s and any B G JF^Z we have 

lim \P(ADB) -P(A)P(B)\ =0 

r— >±oo 

exponentially fast. 

For instance, we can assume that there exists some M > such that P(Af] B) = 
P(A)P(B) for |r| > M. 

Here and below the symbols P and E etc. refer to probabilities and expectations 
etc. with respect to the filtration {J 7 5}- o<s<t<oo- 

In many problems it is natural to assume that the domain D is a thin and long 
channel. This leads to the formulation of the problem as follows. Let D £ = {(x,ez) : 
(x, z) G D}. The parameter e > is small. Consider the diffusion process X t = (Xf , Zf) 
in the domain D £ , which is described by the following system of stochastic differential 
equations: 

dXf = dWl + V(XI, Zf/e)dt + ^(Xf, Zf)Mt , 
dZt = dW? + v 2 (Xf,Zt)dll . 



2 




Fig. 1: A model of the molecular motor. 



Here the scalar field V(x, z) > 0, (x, z) G D characterizes the speed of the trans- 
portation in the x-direction. The vector field v = (^1,^2) on dD e is defined as the 
inward unit normal vector at the corresponding point on dD: f(x, ez) = n(x, z) when 
(x,z) G dD. The process {W t , W 2 ) is a standard 2-dimensional Wiener process. We 
make another assumption here. 

Assumption 5. The process (W/, W 2 ) is independent of the filtration {J-^}-oo<s<t< 
corresponding to the shape of D. 

- — -e 

In other words our process X t is moving in an independent random environment 
characterized by random shape of the domain D. 

The process l\ is the local time of the process X t at dD £ . 

Making a change of variable Zf — >■ Zf/e = Zf in the equation (1), we can equiva- 
lently consider the diffusion process Xf = (Xf, Zf) in the original domain D as follows: 

dXI = dWl + V(Xf, Zf)dt + uf(XI, Zt)di\ , 
dZt = \dWl + vl{Xt,Zf)dl% , (2) 

id 2 id 2 

Here v e is the co- normal vector field corresponding to the operator — ^ 7T~2 + — 7T" 2 : 
u e (x,z) = (vf(x,z),V2(x,z)) = (en\(x, z),n2(x, z)). The process l\ is the local time of 
the process X\ at dD. 

Here and below we use the symbols P w and E w etc. (sometimes with a subscript 
to denote the starting point of the process) to refer to probabilities and expectations 
etc. with respect to the filtration generated by the 2-d Wiener process (Wf,Wf). 

The process X\ has the "fast" and the "slow" components. The "fast" component is 
the process Z\ and the "slow" component is the process Xf. According to the averaging 
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principle we can expect a mixing in the "fast" component before the "slow" component 
Xf changes significantly. We shall describe the limiting slow motion. 

In the next section we will characterize the limiting slow motion, which is a diffusion 
process on a graph. This graph corresponds to the domain D. A sketch of the proof of 
this result is in Section 3. 

An interesting question arising in the applications is to calculate the effective speed 
of the particles. In mathematical language this problem can be formulated as follows. 
Let <7 e ((— oo,a]) be the first time that the process Xf, starting from a point xq = 
(xq, zo) G D, hits D n {x = a > 0}. The limit 

lim hm (7£(( - 00 ' a]) 

a-s-oo elO a 

exists in P x Pj^ ZQ \ -probability and can be viewed as the inverse of the average effective 
speed of transportation of the particle inside D. Using the results in Sections 2 and 3 
we can calculate this limit. This is done in Section 4. (In particular, see Theorem 7.) 

In the last Section 5 we mention briefly problems for multidimensional channels, 
for random channels changing in time, and some other generalizations. 

2 The limiting process 

Let us, for the present and for the next section, work with a fixed shape of D. We 
will allow this shape to be random in Section 4. We shall find the limiting slow motion 
of the diffusion process X\ inside D. 

First of all we need to construct from the domain D a graph T (see Fig.l). For 
xq G K let C{xq) = {(x,z) £ D : x = xq} be the cross-section of the domain D with 
the line {x = xq}. The set C(xq) may have several connected components. We identify 
all points in each connected component and the set thus obtained, equipped with the 
natural topology, is homeomorphic to a graph T. We label the edges of this graph T by 
I\, Ik, ... (there might be infinitely many such edges). 

A point y € T can be characterized by two coordinates: the horizontal coordinate 
x, and the discrete coordinate k being the number of the edge in the graph T to which 
the point y belongs. Let the identification mapping be 2) : D — > Y. We note that the 
second coordinate is not chosen in a unique way: for y being an interior vertex Oj of 
the graph V we can take k to be the number of any of the several edges meeting at the 
vertex O,. 

The distance p(yi,y2) between two points y\ = (x±,k) and yi = (x2,k) belonging 
to the same edge of the graph T is defined as p(yi,y2) = \x\ — X2\; for y\,y2 G T 
belonging to different edges of the graph it is defined as the geodesic distance p(yi,y2) = 
min(p(yi, Oj x ) + p(Oj 1 , Oj 2 ) + ... + p(Oj v 2/2)), where the minimum is take over all chains 
yi f-> Oj 1 -H- Oj 2 f> ... O Oj l -H- y2 of vertices Oj. connecting the points y± and 2/2 ■ 
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For an edge I k = {(x,k) : A k < x < B k } we consider the "tube" U k = 2) _1 (ijfc) n 
{A k < x < B k } in D. The "tube" Uk can be characterized by the interval x £ [A k ,B k \ 
and the "height functions" h k (x): Uk = {(x,z) : Ak < x < Bk,h^(x) < z < h^(x)}. 
For x £ [Ak, Bk], we denote the set Ck(x) to be the connected component of C(x) that 
corresponds to the "tube" Uk'- Ck(x) = {x} x [h k ,h^]. Let lk{x) = h~^(x) — h^(x) > 
Ik > for all i£i We notice that each h^{x),lk{x), etc. is smooth. 

The vertices Oj correspond to the connected components containing points {x, z) E 
dD with 712(2,2:) = 0. There are two types of vertices: the interior vertices (in Fig.l 
they are 0\, O3, O5, O7) are the intersection of three edges; the exterior vertices (in Fig.l 
they are O2, 04, Oq, 0%) are the endpoints of only one edge. 

Using the ideas in [H] with a little modification we can establish the weak conver- 
gence of the process Yf = %)(Xf ) (which is not Markov in general) as e 4, in the space 
C[ 0i T](r) to a certain Markov process Y t on F. 

The process Yt is a diffusion process on V with a generator A and the domain of 
definition D{A). 

For each edge Ik we define an operator 

— 1 d ( , du\ — du 

L k u(x) = — lk{x)— + V k {x)— , A k < x < B k . 



2l k (x) dx \ dx J dx 
Here 

V k {x) = y-r -r \ V(x,z)dz 

h(x) J h - {x) 

is the average of the velocity field V(x,z) on the connected component Ck(x), with 
respect to Lebesgue measure in z-direction. 

For simplicity of presentation we will assume throughout this paper the following. 

Assumption 6. The function Vk{x) = ft > is a constant. 

The case of non-constant Vk(x) can be treated in a similar way. The only difference 
is that the calculations are a little bit more bulky. 
We also let 

The operator L° can be represented as a generalized second order differential op- 
erator (see [T]) 

lP k u(x) = D rk D g J(x) , 

Q\X -j - OJ — Cjix) 

where, for an increasing function h, the derivative is defined by D^gix) = lim — , 

(5J.0 h{x + 0) — h(x) 

and 



Qk{x) 



dx f 
— - , r k (x) = 2 1 l k (x)dx . 
h{x) J 
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The operator A is acting on functions / on the graph T: for y = (x, k) being an 
interior point of the edge Ik we take Af(y) = Lkf(x, k). 

The domain of definition D (A) of the operator A consists of such functions / 
satisfying the following properties. 

• The function / must be a continuous function that is twice continuously differ- 
entiable in x in the interior part of every edge Ik', 

• There exist finite limits lim Af(y) (which are taken as the value of the function 

Af at the point Oi); 

• There exist finite one-sided limits lim D qk f(x, k) along every edge ending at 
Oi = (xi, k) and they satisfy the gluing conditions 

Ni 

where the sign "+" is taken if the values of x for points (x, kj) G are > X\ and "— " 
otherwise. Here N% = 1 (when Oi is an exterior vertex) or 3 (when Oi is an interior 
vertex) . 

For an exterior vertex Oi = (xi , k) with only one edge Ik attached to it the condition 
(3) is just lim D qk f(x, k) = 0: the boundary condition of instantaneous reflection at the 
end Oi. (We notice that we are in dimension 2, in dimension 3 or higher these endpoints 
do not need a boundary condition, they are just inaccessible.) For an interior vertex the 
gluing condition (3) can be written with the derivatives — instead of D qk . For k being 
one of the kj we define ctik = lim lk(%) (for each edge Ik the limit is a one-sided one). 
Then the condition (3) can be written as 

^-^ J x-tXi dx 

i=i 

It can be shown as in |5j Section 2] that the process Yj exists as a continuous strong 
Markov process on T. 

For every e > 0, every x = (x, z) £ D and every T G (0, oo) let us consider the 
distribution [i e x of the trajectory Y t £ = %)(Xf ) starting from a point Xq = x in the space 
C[o ) T](r) of continuous functions on the interval [0, T] with values in T: the probability 
measure defined for every Borel subset B C C[ 0r ](r) as n%{B) = Pxg=a;(^» £ e B). 
Similarly, for every y G V and T > let fj,y be the distribution of the process Yj in the 
same space: \r v (B) = P^{Y. G B}. The following theorem is our main tool for the 
analysis. 

Theorem 1. For every x G D and every T > the distribution \i e x converges 
weakly to to $qi x \ as e \,0. 
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In other words we have 



for every bounded continuous functional F on the space C[o ; T](r)- 



3 Sketch of the proof of the convergence 

We shall now briefly give a sketch of the proof of Theorem 1 announced in the 
previous section. 

The averaging within each edge Ik is a routine adaptation of the arguments of [5j 
Section 3]. Within one edge Ik, the motion of the component Xf is given by the integral 
form of the stochastic differential equation 

X £ t = W} + / V(XI, Zt)ds + / uf(XI, Zl)M% , 
Jo Jo 

and the one for the limiting motion Xt looks like 



From the above two formulas we see that in order to prove the convergence of X 
to Xt as e \. in the interval < i < T < oo we just need the estimates of 



and 



{I) 



(II) 



AO z _ I [* l 'k( X !) . 

* H i k (xi) ds 



Jo 



f v(xi,zi)ds- fvkixDdt 

Jo Jo 



The estimate of (I) ("averaging with respect to local time") is exactly the same as 
that of [5l Section 3]. For the estimate of (II) we can introduce an auxiliary function 
cik(x,z) satisfying the problem 



1 d 2 a k 

2 8z 2 



(x,z) = V(x,z) - V k (x) , 



. ^,K(x)) = ^(x,ht(x))=0. 
The solvability of this equation is guaranteed by the fact that 



/ 

Jh 



{V(x,z)-V k (x))dz = 



h k (x) 

(this is the key point in averaging). The solution is bounded with bounded derivatives. 
Applying the generalized Ito's formula (see Section 3, equation (3.1)]) to the function 
dfc we see that 
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a k (Xf,Zt)-a k (X' a ,Zl) 

:' + 1 £ to* + jC 

Multiplying both sides by e and taking into account the problem that satisfies 
it is immediate to get an estimate of (//). These justify the averaging within one edge 
h- 

The gluing conditions can be obtained using the results of [5] and the Girsanov 
formula. To this end one can introduce an auxiliary process X t = (Xf,Zf) in D via 
the following stochastic differential equation 

dxi = dwl + uf(Xf,Z!)dit , 

dZf = - £ dWf + ul{Xt,ZI)dlt . 
Here i\ is the local time for the process X t at dD. 

This is exactly [51 formula (1.4)]. The limiting process within an edge is governed 
by L^. The gluing condition is just the gluing condition in (3). On the other hand, The 
measure m e corresponding to the process X\ is related to the measure rh £ corresponding 
to the process X t in C\q^\(D) via the Girsanov formula 



:{X\) =e W S [ £v(X £ s )dW} - ±£v 2 (Xl)ds\ 



dm £ ~e, f /" T „,~e s „„i 1 

dm 6 

The Girsanov formula and a standard argument (see [3] and [31 Appendix A. 2]) 
guarantee that the gluing conditions remain the same for the process X\ and thus we 
have proved Theorem 1. 



4 Analysis of the limiting process 

The next goal is to analyze some interesting quantitative characteristics of the 
one-dimensional diffusion process Yt on T. As we have pointed out in Section 1, the 
interesting question coming from applications is to quantify the effective speed of the 
motion of the particle. Our calculation in this section will always be performed by first 
fix a shape of D and then let the shape of D be random. 

4.1 The case when there is only one edge of the graph T 

The simplest case is that when there are no "wings" and also that the graph T 
consists of only one edge I\. Let the corresponding h(x), V\{x) = (3 > 0, etc. be 
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defined. Recalling our assumptions in Section 1, we see that in this case the functions 
h (x) do not have discontinuities and lo(x) = h(x) is smooth and positive, uniformly 
bounded from above and below. Consider the interval (—00, a] for some a > 0. Let the 
process Yt start from x = 0. Let r((— 00, a]) be the first exit time of Yt from (—00, a]. 
The random variable lo(x) = h(x) is distributed according to our stationary and mixing 
assumptions. We have the following. 
Theorem 2. We have 



a— >oo a 



K(t) exp(-2(3t)dt 



where the function K(t) = E- 



l (s + t) 

Proof. We see that u(x) = Ej^r((— 00, a]) is the solution of the problem 



1 d f du\ ^ ^du 
2li(x) dx \ dx J dx 



■1 , 



u{— 00) = n(a) = . 

The above problem can be solved explicitly. We shall first expand the equation as 



Now we introduce 



dx 2 



+ 



l'i(x) ^ du 
h(x) J dx 



-2 . 



<p(x) 
= exp 
_ h(x) 



x (I'M 
h(y) 



+ 2/3 )dy 



du , 
dx 



du s 

Zi(o) exp(2/3x)-(,) 



Using the equation that u(x) satisfies it is not hard to check that 



Integration gives 



Thus 



<p(x) 



h(0) 

h(y) 



exp (2/3x) 



h(0) 



exp (2py) dy + y>(0) 



^ = -t^t r^i(y)exp(-2/3(x- ?/ ))d 2 / + v3 (0)^exp(-2/3x) . 
Taking into account that u(a) = we see that 
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u(x) 



dy 

a h(y) \Jo 



h (t) exp (-2/3(y - t)) dt\ + (^(0) J ^ exp (-2(3y) dy . 



Thus 



»< > = 2 frrr, 

Jo h{y) Jo 



h(t) exp (-20G/ - t)) - C(a) 



1 



o z i(y) 



exp (-2/3y) . 



Here the constant C(a) can be determined from the fact that lim u{—b) = 0. This 

fe— S-oo 



gives 



C(a) 



lim 

b— >oo 



lim 

ft— »oo 



a dy 



1 



My) 



exp (-2/3y) dy 



«i(-6)7o 



Zi(t) exp (2/3(6 + t))dt 



1 



exp (2/36) 



-2 / h(t) exp (2/3t) dt . 

J — oo 



We see from above that |C(o)| < C < oo, where the constant C is independent of 
a. This, combined with the fact that 



1 



■ exp (-2/%) dy 



'o *i(y) 

is uniformly bounded in a, show that the limit is equal to 



a— >oo a 

r 2 f a dy 
= urn — 



h(t) exp(-2/3(y - i))dt 



a Jo h(y) \Jo 
2 f a dv I f y 
lim - / / h(y-t)exp(-2pt)dt 

9 /-a /-a-s i f \ 

lim - ds p ' - exp(-2/3t)dt . 

a^ooaj / /i(s + t) 



Let the random variable L(s) 



, exp(-2/3t)dt. We have \L(s)\ < L is 
/o «i(s + *) 

uniformly bounded. We fix an arbitrary < \i < 1 and we have, for any k > and any 
< s < fxa, there exist ag = ao(a, K, //) > such that for any a > «o we have 



h(s) 



o Ms + t) 



exp(-2(3t)dt - L(s) 



< kL . 
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Thus 



o ffia t*a—s 

- ds 

a Jo Jo h{s + t) 



U ' exp(-20t)dt - - / L(s)ds 
a Jo 



< 2fiK,L 



On the other hand, by our mixing and stationarity assumptions we see that there 
exist Qi = a\(n, /i) > such that for any a > a\ we have 



2 f^ a 
a 

Therefore when a > max(ao,ai) we have 



/ L(s)ds - 2fiEL(s) 
Jo 



< K . 



h(s + t) 
On the other hand we have 



2 - r d3 

a Jo Jo 



exp(-2pt)dt - 2fiEL(s] 



< 2/xkL + k . 



a pa~s 

ds 



fia 



h(s) 
h{s + t) 



exp(-2(3t)dt 



< 2L(1 - fj) . 



Thus when a > max(ao, a±) we have 



o pa ra—s 

t ds 
« Jo Jo 



h(s) 
h(s + t 



exp(-2(3t)dt - 2BL{s] 



< 2(inL + k + 4L(1 - fj,) . 



Since we can take [i arbitrarily close to 1 we see from the above estimate that we 



have 



a^oo a Jo J 



exp(-2f3t)dt = 2BL(s) . 



h(s + t) 

Taking into account that lo(x) = l\{x) in this case we conclude with the statement of 
this theorem. □ 



4.2 The motion inside D 

Now we consider in more detail the case when the graph T does not have any branch- 
ing but it has many edges I±, Ik, ••• that form a straight line. Let Oi, Ok, ■■■ etc. be 
the corresponding vertices. Let the corresponding cross-section width be l±(x), ...,lk(x), ... 
etc.. (See Fig. 2.) In this case the functions h (x) have jumps. We can introduce a func- 
tion Iq(x) = lk(x) for Ak < x < -Bfc. The function lo(x) has jumps at Oj's and at the 
jumps of lo(x) we connect the pieces of the boundary via vertical straight lines. In this 
way we form the domain Dq as we introduced in Section 1. We can find a family of 
smooth functions h^' S (x) such that h,Q ,S (x) converge as 5 1 uniformly on compact 
subsets of R to h^(x). Thus h^' S {x) — ,S (x) = Iq(x) converge as S 1 uniformly on 
compact subsets of R to the function Iq(x). Also, we can choose such ^^(x) that the 
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domain Dq = {(x, z) : h^ ,S < z < h Q ,S } satisfies our assumptions. Consider the process 
Xf in the domain Dq defined as in (2). Let Y t £ '° = Z)(X £ t '°). For fixed 5 > 0, the 
graph r corresponding to Y^' S has only one edge I\. Using the results in Section 2 we 
see that ai 
generator 



see that as e J, the processes Yf' S converge weakly to Y t s . The process Y t s on V has a 



L ^) = ^7-T3Z fez )+V 1 (x)- 



Here 



tt5 



2l s Q (x) dx \ 0K ' dx J 1 dx 



^iW = ^7T t V(x,z)dz 



and we have V 1 (x) — > Vk{x) = (3 > as 5 — > uniformly on compact subsets of x E R. 

The above operator can be written in the form of a D v sD u s operator in the sense 
of PQ. We can explicitly calculate the functions u s and v s as follows: 



uHx) =£wM- 2 L ,vt(t)dt ) dv 



v°(x) = 2 J ^(y)exp ^2 J Vl{t)dt\ dy . 

Consider a Markov process Y t D ° on the graph T. The process Y t D ° is governed by 
the generalized second order differential operator L = D V D U (in the sense of pQ) where 



u(x) = J exp ^-2 jf " V k (t)dtj dy 



u 



v(x) = 2 J l (y)exp \2 J V k (t)dt\ dy . 

The domain of definition of the operator L consists of those functions f(x) that are 
continuous and bounded, are twice continuously differentiable and at those points 
where lo(x) have jumps it satisfies a gluing condition j^rrfLiQi) = J7r\f+{0i). 

Here lo t ±(Oi) are the corresponding values of the one sided- limits of Iq(x) to the left and 
to the right of Oi, and /4(Oi) the corresponding left and right derivatives. The function 
Lf{x) is continuous on T. 

It follows from the classical result of [6] that we have the following. 

Theorem 3. As 5 1 the processes Yj* converge weakly to the process Y^° on 
C[o,T](r). 

Let t <5 ((— oo,a]) be the first time when the process Y t s exits from (— oo,a]. Let 
T D °((—oo,a\) be the first time when the process Y t D ° exits from (— oo, a]. Weak con- 



vergence of processes Y t s as 5 I to Y t D ° implies that we have limEQ / r <5 ((— oo, a]) 

(54,0 

E^((-oo,a]). 
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We recall the first differential equation we used in the proof of Theorem 2. We 
plug in l\(x) = Iq(x) in that equation and we see that the corresponding solution is just 
the solution we get there with l\{x) replaced by lf(x). However Eq / t <5 ((— oo, a]) is the 
solution of the same problem we used in the proof of Theorem 2 with h(x) replaced by 
Iq(x) and (3 replaced by v\{x). Since v\{x) — > j3 as 5 1 we see that we have 



E?V*»((-oo,a]) = limE *V((-oc,a]) 

o\.0 



hm (2 P ^- r l 5 (t) exp (-2/3(2/ - t)) dt + 2 f° l s (t) exp (2(3t) dt f J- exp (-2/3y) 

^° V Jo w) Jo J-00 Jo l {y) 

lo(t) exp (-2/3(1/ - *)) dt + 2 / fo(t) exp (2pt) dt [ exp (-2(3y) dy . 

J -00 Jo hw 



dy 
lo(y) 



Following our stationarity and mixing assumptions, after deleting the "wings", the 
remaining channel still satisfies the stationarity and mixing assumptions. Therefore by 
the same calculation as in the proof of Theorem 2, and using the above formula, we see 
that we have the following. 

Theorem 4. For the process Y t defined as above we have 

lin E^((-oc, a| ) =2 ^ 

a->oo a Jq 

where Kit) = and we allow jumps of the function Iq(x). 

h{s + 1) 

On the other hand, using the same argument of [5] one can show that the process 
Y t D ° can be viewed as the limiting slow motion as e J, of the part of the process X £ t 
within the domain Dq. To be precise, consider the domain D introduced in Section 1 

and the corresponding process X £ t in D. Let <ft = l l(-^s £ Do)ds be an additive 

J 

functional. (It is called the proper time of the domain Do, see [8].) We introduce the 
time Pt inverse to cpt and continuous on the right. Let Y^' D ° = %)(X% t ). Then we can 
use the same arguments of [5] to prove the weak convergence as e ! of the processes 
y/' D ° to the process Y t D °. 

Consider the process Xf moving in the domain D as in Section 1. Let Iq(x) be 
the cross-section width corresponding to the domain Dq. As before we see that it could 
have jumps. Let Yt be the limiting slow motion as e I of Xf, defined as in Section 
2. Let t((— 00, a]) be the first time the process Yt exits from (—00, a]. Since Yt is the 
limiting slow motion of X\ and Y t D ° is the limiting slow motion of the part of X\ inside 
Dq, we see from the above discussions that we have 

/■r((-oo,o]) 

T D °((-oc,a})= l(fQ-\Y t ) C D )dt . 

Jo 

By Theorem 4 and the above relation we see that we have the following. 



13 




Fig. 2: The case when Iq(x) has jumps. 



Corollary 1. 

/•r((-oo,a]) 

E ^/ lty-^Yt) C D )dt 

lim ^ = 2 / K{t)exp(-2f3t)dt 

a Jo 



a— >oo 



w/iere -fT(t) = E /o(s 



l (s + t)' 

This fact will be used in the next subsection. 



4.3 The general case 

In the general case when there is branching the domain D consists of a domain Dq 
that has a cross-section width Iq{x) which allows occasional jumps. The "wings" Dj 
(j > 1) are then attached at the jumps of the domain Dq. 

In this case our graph V consists of two types of edges. The first type of edges 
correspond to the domain Dq as we discussed in the previous section. The second type 
of edges correspond to the "wings" attached to Dq. 

In order to calculate the effective speed of transportation we shall first calculate 
the expected time that the process spends at one fixed "wing" . As a first step we do not 
consider the random shape but perform the calculation for a fixed shape. Also, we shall 
first consider the simplest case that T has only three edges: I\ = (— oo,0], h = [0,oo) 
and ^3 is an edge with one endpoint 0\ = and another endpoint O2 = (^,3). In this 
case h(x),h{x) and h(x) are smooth functions. 

We construct the process Yt as in Section 2 corresponding to the above graph T. 
Consider the interval (—00, a] for some a > 0. Let the process Yt start from the point 
0\ = 0. Let tj 3 ((— 00, a]) be the time that the process Y t spends at the edge ^3 before 
its first exit time from (—00, a]. We have the following. 

Lemma 1. 
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E^r /3 ((-oo,a]) = 2sign(r) f l 3 (t)exp(2pt)dt [ -L- exp 

Jo Jo l 2{y) 



Proof. Let for example r > 0. The function u(x,k) = E|^ fc ^r/ 3 ((— oo, a]) is the 



solution of the problem 
f 1 d 



i l k(x)p- J + = , for (x, fc) G 4 and fc = 1, 2 , 
2lk(x) ax \ dx J dx 

1 d / ' , .du\ „du , 

it(— oo, 1) = u(a, 2) = , 
diti . . 

Zi(0) lim ^(x,l)=Z 2 (0) lim ^(x, 2) + fe(0) lim ^(x,3) . 
i->Oi dx x^Oi dx x-^Oi ax 

The function u(x,k) is continuous at point 0\. 

Similarly as in Section 4.1 there are solutions u\[x) = u(x, 1), 7/2(x) = u(x,2) and 
7/3 (x) = u(x, 3) corresponding to the edges Ji, J2 and I3. They are defined as follows. 
We have 



X 



u 1 {x) = C 1 h(0) 

U2 (x) = cy 2 (o) 



1 



b h(t) 
x 1 



a h{y) 



exp(-2(3t)dt + ui(-b) for x G (-00, 0] , 
exp(—2/3y)dy for x G [0, +00) , 



1 



exp(-2/3y)dy + 7/3(7-) 



h(y) 

for x G [0, r] 



We also know that lim 7/1 (— 6) = 0. There are 4 undetermined constants: C\,C2,C 3 
and 7/3(7*). We can uniquely determine them by the 4 relations 

7/i(0) = 7/2(0) = 7/3(0) , 

^3/ \ n 
ax" (r) = ' 



and 



-M0)^(0) + M0)^(0) + 1 3(0)^(0) = 0. 

We shall solve the above system. The relations we need are as follows 
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1 

-b h{y) 



-C 2 Z 2 (0) I -^-exp(-2py)dy 

o V ' o 
-2 J ^n% 3 (t)exp(2p(t-y))dt\+C 3 l 3 (0)J j^-r exp(-2/3y)dy + u 3 (r) 

rr 

-2 / Z 3 (i) exp(2/3(t - r))dt + C 3 Z 3 (0) exp(-2/3r) = , 
-C 1 l 1 (0) + C 2 h(0) + C 3 l 3 (0) =0 . 

Thus 

C 3 Z 3 (0) = 2 / l 3 (t)exp{2pt)dt . 
Jo 

On the other hand, we have 

dhtp) = C 2 / 2 (0) + C 3 / 3 (0) . 

So we get 



(C 2 / 2 (0) + C 3 / 3 (0)) ' 



/ -l_exp(-2/3y)dy + «i(-6) = -C 2 Z 2 (0) / exp{-2/3y)dy . 

J-b h{y) Jo h{y) 



This gives us lim (C 2 Z 2 (0) + C 3 Z 3 (0)) = 0. Thus 

b— >oo 

u 2 (0) = 2 / l 3 (t)exp(2pt)dt [ -J—exp(-2/3y)d2/ . 
Jo Jo «2(y) 

In the case when r < the gluing condition becomes 

, / s du , , , , . , du , , , , . , dit . 
Zi(0) lim — x,l) + Z 3 (0) lim — (x, 3) = Z 2 (0) lim — x,2 . 
x->Oi ax nOi ax i^Oi ax 

Thus CiZi(O) + C 3 / 3 (0) = C 2 Z 2 (0). A similar calculation shows that 

E y r /3 ((-oo,a]) = 2 [ l 3 {t)exp(2(3t)dt [ -L- exp(-2py)dy . 

Jr JO [ 2{y) 

Thus in general we have 

E^r j3 ((-oo,a]) = 2sign(r) / l 3 (t) exp(2/3t)dt / e xp(-2(3y)dy . 

Jo Jo h(.y) 

□ 

The above lemma can help us to deal with a more general case. Let q ^ and 
q £ (—00,0]. We assume that the graph T still consists of 3 edges I\ = (— 00, q], 
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li = [q,oo) and I3. In this case the edge I3 has one endpoint 0\ lying on I\ U I2 but 
with x-coordinate q, and another endpoint O2 = (q + r, 3). Let the process Y t again 
start from the point 0. (In this case the point is lying on either I\ or I2 but may not 
be on their intersection.) Let t/ 3 ((— 00, a]) be the time that the process Yt spends at 73 
before it exits from (—00, a]. We have the following. 
Corollary 2. If q > then 

E^r/ 3 ((-oo,a]) = 2sign(r) J* + l 3 (t) exp(2/3(t - q))dt J j-^ exp(-2/% - q))dy . 

Ifq<0 then 

E^r /3 ((-oo,a]) = 2sign(r) T + l 3 (t) exp(2/3(i - g))di / — L e xp(-2/% - q))dy . 

Jq Jo h{y) 



Proof. The above results can be easily seen from the strong Markov property of the 
process Yt and Lemma 1. To be more precise, if q > 0, then by strong Markov property 
Eq / t/ 3 ((— 00, a]) is just equal to E^ / r/ 3 ((— 00, a]), which can be calculated using Lemma 

1 and a shift. If q < we need a little bit more argument. In this case we can consider the 
process Yt starting from x = q and its first time of exiting from (—00, 0]. Let r/ 3 ((— 00, 0]) 
be the time that Yt spends at I3 before it exits from (—00, 0]. Then E^ / r/ 3 ((— 00, 0]) can 
be calculated using Lemma 1 and a shift. On the other hand, E^ / t/ 3 ((— 00, a]) can also 
be calculated using Lemma 1 and a shift. From the strong Markov property of Yt we 
see that Eq / t/ 3 ((— 00, a]) = E^ / r/ 3 ((— 00, a]) — E^ / r/ 3 ((— 00, 0]), which gives the formula 
we need. □ 

Following a similar approximation argument as we did in Section 4.2, we can con- 
sider the case when the graph T consists of many edges I\,l2, ■■■ etc. that are of the 
first type. They correspond to the domain Dq. We allow jumps of the function Iq(x). 
Then we attach only one "wing" D\ to Dq and the domain D\ corresponds to an edge 
J in the graph T. Let / w ing(^) ; r etc. be the quantities corresponding to the "wing" 
D\. Let the edge J have two endpoints: Of lying on US^lfc and with x-coordinate q; 

02 lying on the other endpoint of J with x-coordinate q + r. After an approximation 
and averaging we get a process Y t D ° uDl on T. Let positive a > q. Let tj((— 00, a]) be 
the time that the process Y t D ° uDl spends in the edge J before it exits from (—00, a]. 
Lemma 1 and Corollary 2 and the same approximation argument as in Section 4.2 give 
us the following corollary. 

Corollary 3. If q > then 

E£Vj((-oo, a}) = 2sign(r) J" + l whlg {t) exp(2/3(t - q))dt J ^-y exp(-2/3(y - q))dy . 
If q < then 

E^rj((-oo, a]) = 2sign(r) T + Z wing (i) exp(2/3(t - q))dt [ -I- exp(-2/3(y - q))dy . 

Jq Jo h(y) 
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If q = then 

Eo / rj((-oo, a]) = 2sign(r) / l wing (t) exp(2(3t)dt / — - exp(-2/3y)dy . 

Jo Jo l o{y) 

Finally we come to the original problem in Section 1. We consider the case when 
there are many random "wings" attached to -Do- Let the process corresponding to the 
domain D be X.% . Let the limiting slow motion be Yj. We introduce the corresponding 
quantities l w - m g(x), r etc. Let Iq(x) = ^q(x) — hg(x) be the cross section width that 
corresponds to the domain Do- We notice that it has occasional jumps. Let there be n 
wings in the interval x G [0, 1]. We assume that we have the following. 

Assumption 7. The random variable n is a bounded random variable: \n\ < < 

oo. 

For a > we define the random variable 



2sign(r) / l w - mg {t) exp(2/3(t - q))dt / — - exp(-2/% - q))dy when q > 

Jq Jq h{y) 



M(q,r,a) 



2sign(r) j ^wing (t)exp(2/3(t - q))dt 



2sign(r) / l wing (t) exp(2/3t)dt 



1 



lo(y) 



exp(— 2(3(y — q))dy when q < 



h(y) 



exp(—20y)dy when q = 



We assume that, for some constants < A$ < A\ < oo, we have the following. 
Assumption 8. Aq < l w i ng (x) < A\ and |r| < A\. 
Thus we see that for some constant A > we have 



^wing(i) exp(2(3t)dt 



< A < oo . 



We also see that 



r°° i 

\M(q,r,a)\ <2A - exp(-2^y)dy < M < oo 
Jo 'o 

for some constant M > 0. 

All these random quantities are distributed according to our stationarity and mixing 
assumptions. Let r((— oo,a]) be the first time the process Y t exits from (— oo,a]. 

By our assumption on stationarity and mixing, using Corollary 2, we see that we 
have the following. 

Lemma 2. We have 



&0 



lim ■ 

a— >oo 



t(— oo,a] 



m-\Y t ) <t D )dt 



2E?iEsign(r) / l wing (t) exp(2/3t)dt 



1 



o l o{y) 



exp(-2(3y)dy 



in probability. 
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Proof. Let the "wings" located to the left of the point have x-coordinate > 
q-i > g_2 > •••• Let the "wings" located to the right of the point (including possibly 
the point 0) and not exceeding x = a have x-coordinate < q\ < ... < q n ( a ) < a. We 
see that we have 

pt(-oo,o,] 00 n ( a ) 

/ l(r 1 (r 1 )fO ) ( lt = VM(9_ l ,r_ t ,fl)lVM( ft ,r l ,fl) . 

Jo k=i k=i 

For q < we have 

M(q,r,a) < exp(2/3q)M(0,r,a) < M exp(2/3g) . 

Thus 



oo 

Y,M(q-k,r- k ,a) < M^exp(2/3g_ fcJ 

k=l k=l 

Taking into account Assumption 7 we see that we have 

oo oo 

exp(2/3g_fc) < np exp(— 2f3k) < oo 

k=l k=0 

Therefore 



/■r(— oo,o] n(a) n(a) 

Kj HV\Y t ) £ D )dt £M(« fc ,r fc ,a) "£M(q k ,r k ,a) 
lim = lim — = lim 



a— >oo a a— >oo a a— too n(d) 

On the other hand, we have 

n(a) 

,r k ,a) 

k=l 

= z2 2sign(r fc ) / l W mg(t) exp(2/3(i - q k ))dt I —— exp(-2/3(y - q k ))dy 

k= i Jqk -jQk l o\y) 

n(a) 

= 2Y,W(q k ,r k )C(q k ,a) . 
k=l 

Here 

rik+rk 

W(q k , r k ) = sign(r fc ) / Z w ing(*) exp(2^(t - q k ))dt 

•><lk 

and 

f a 1 

C{q k , a) = / — - exp(-2/3(y - q k ))dy . 



<ik 



. h{y) 
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Thus we can write 

n(a) n(a) n(a) 

^2 M(q k , r k , a) =2^2 W(q k , r k )C{q k , a+q k )-2^2 W (Vk, r k ){C(q k , a+q k )-C(q k , a)) . 

k=l k=l k=l 

By the remark after Assumption 8 we see that 

n(a) 

2 ^2 W (lk, r k )(C{q k , a + q k )- C{q k , a)) 
k=l 

n(a) 

< 2A J2(C(q k ,a + q k ) - C{q k ,a)) 
k=l 

2A ^> r a+q k 



2A v- 1 ^ r a+ ik 
^ ~T Yl / e M-W(y - Qk))dy 

[ o 7~, J a 



k=l 

^[exp(-2p> - q k )) - exp(-2/3a)] 

k=l 



A ^ A 
= — V exp(-2/3(a - q k )) - -^-n(a) exp(-2/3a) . 

Thus by our Assumption 7 again we see that 

n(a) 

2 X) W(g/ c ,^fc)(C , (g fc ,a + g , fc) -C(q k ,a)) 
lim — = . 



a— ^oo 



n(a) 



Therefore using the weak Law of Large Numbers for triangular arrays and taking 
into account our assumptions on mixing and stationarity we see that 

n(a) 

X M (lk,r k ,a) 
k=l 



lim 

a-^oo n[a) 



n(a) 

E W(q k ,r k )C(q k ,a + q k ) 
k=l 



= 2 lim 

a->oo n(a) 

= 2Esign(r)/ / wing (t)exp(2/3t)dt / — — exp(-2/3y)dy 

Jo Jo l o(y) 

in probability. 

On the other hand we have lim = En almost surely. Thus we can conclude 

a— >oo a 

with the final result. □ 

Adding the two equations in Corollary 1 and Lemma 2 we have the following. 
Theorem 5. We have 
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Hm EgV((-oo,a]) 







pr poo 
Kit) exp(-2f3t)dt + 2EnEsign(r) / l win Jt) exp(2f3t)dt / — - exp(-2(3y)dy 

Jo Jo low 

in probability. Here Kit) = E ^ ^ . 

l (s + 1) 

Let cr e ((— oo, a]) be the first time that the process Xf, starting from a point Xq = 
(xq,zq) inside D, hits the level curve C(a). Since xq is finite we see that Theorems 1 
and 5 lead to the following theorem. 

Theorem 6. We have 

y r E^^((-oo,q]) 
lim lim — - 

a-^oo elO a 

/•oo rr roc i 

= 2 / Kit) exp(-2/3t)dt + 2EnEsign(r) / l win Jt) exp(2f3t)dt / — - exp(-2(3y)dy 
Jo Jo Jo h{y) 

in probability. Here Kit) = E — \° . 

h(s + t) 

The above theorem helps us to conclude in our original problem. We can divide 
the domain D n {{x,z) : x G (— oo,a]} into consecutive pieces alternatively of x-length 
b and y/b for some b > 0. The total time spent by the process X\ before it exits from 
x = a is the sum of those times spent in domains of x-length b and those times spent in 
domains of x-length \fb. As we are taking a — > oo we can also let b — > oo and the average 
time spent in domains of x-length y/b will not contribute. On the other hand, since the 
process X\ has a deterministic positive drift in the x-direction, we see that as b — > oo 
the motion inside different domains of x-length b will be asymptotically independent. 
(The motion against the flow is a large deviation effect.) Thus a e ((—oo,a\) will be 
asymptotically distributed as the sum of some independent random times. This leads 

to the relation lim lim — — — — ■ — — = lim lim — — — — . Thus from Theorem 6 

a-Kx> eiO a a->oo e i0 a 

we see that we have the following theorem. 
Theorem 7. We have 

lim lim CT£(( ~°°' a]) 

poo pT poo 

= 2 1 K(t)exp(-2f3t)dt + 2EnEsign(r) / l wing {t) exp(2/3t)dt / — - exp(-2/3y)dy 

Jo Jo Jo io{y) 

in probability. Here Kit) = E ° ^ ^ . 

l {s + 1) 
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5 Remarks and generalizations 



1. The results of previous sections can be extended to the case of a multidimensional 
channel G = {(x, z) G x } where G x C R n_1 , are (n — l)-dimensional 
domains assumed to be bounded and consisting of a finite number of connected com- 
ponents; each G x contains the origin. The boundary dG of G C R n is assumed to be 
smooth, except, maybe, a number of (n — 2)-dimensional manifolds (like the discontinu- 
ity points in the 2-dimensional case). Let T be the graph homeomorphic (in the natural 
topology) to the set of connected components of the sets {z E R™" 1 : (x, z) E G} for 
various x E R 1 . Let the functions k(x) be defined as (n — l)-dimensional volumes of 
corresponding connected components. 

d 1 

Consider the process Xf governed by the operator b— — | — A x z inside the domain 

ox 2 

G £ = {(x,z) : x E R\z E R n_1 , (x,z/e) E G}. Let a £ a = min{t : X% E {(x,z) E 
G U dG,x = a}}. Then, under mild additional conditions the limit 

lim lim — 

a->oo eiO a 

exists in probability P x P|^ z ), (x, z) E G £ and is given by Theorem 7. 

2. Let vt be a continuous time Markov chain with 2 states 1 and 2. Let functions 
hf(x), x E R 1 , i = 1,2, be piecewise smooth, and h~(x) < < hf(x), x E R 1 . Put 
Df = {(x,z) : x E R 1 ,— h^(x) < z/e < hf(x)} and D £ (t) = Dlny Define the process 

d 1 

(X £ , Z £ ) in D £ (t) as the process governed by the operator b— — h — A XjZ inside D £ (t) 
with the normal reflection on the boundary at the times when Uf is continuous (then vt 
is constant). Let Zf jump to at times when vt has jumps. (Actually, we need this 
condition to define the process (X^, Zf) in a unique way; it is not important since we are 
interested in the limit as e 4- 0.) Then one can prove that the slow component Xf of the 
process converges as e I to the process described by the equation X t = Wt + ^H*' ^t) 

where b(t,x) = b + — (lnZ^(z)), l Vt (x) = (x) — h~ (x). It is known that for certain 

_ ax 

drift terms b(t,x), the process Xt demonstrates the so called ratchet effect: If Ut is 
identically equal to 1 or 2, the process tends to +oo as t — > oo, but if ut is the Markov 
chain (independent of Wj) the process X t tends to — oo. 

One can conclude from our considerations, that the ratchet effect can be caused by 
random and independent of the basic Brownian motion changes of the geometry of the 
domain. 

3. We assumed that the "wings" have a simple structure - each of them corresponds 
to just one edge of the graph (Fig.l). One can consider the case of more complicated 
wings, like, for instance, at the vertex 0\ in Fig. 3. One can also include in the consid- 
eration the case when the "obstacles" in the channel are such that the corresponding 
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Fig. 3: A more general graph. 

graph has loops like that in Fig. 3. 

If the channel D e is not "uniformly narrow" but has points on axis x such that in 

the (^-neighborhoods, 5 = 5(e) « 1 of those points the channel has the "diameter" of 

order //(e) » e, the limiting process on the graph can have delays or even traps. This 

a 6 

will lead to different behavior of — . 

a 
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